High three dimensional thermoelectric performance from low dimensional bands 
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Reduced dimensionality has long been regarded as an important strategy for increasing thermoelectric per- 
formance, for example in superlattices and other engineered structures. Here we point out and illustrate by 
examples that three dimensional bulk materials can be made to behave as if they were two dimensional from the 
point of view of thermoelectric performance. Implications for the discovery of new practical thermoelectrics are 
discussed. 



INTRODUCTION 

Thermoelectric performance is quantified by the figure of 
merit, ZT = oS 2 T /k, where a is the electrical conductivity, 
K is the thermal conductivity, S is the thermopower (Seebeck 
coefficient) and T is the absolute temperature. Ill |2[] There 
is no known thermodynamic or other fundamental limitation 
on ZT, but finding high ZT materials is very challenging and 
only a handful of materials with ZT significantly higher than 
unity are known. The difficulty is that finding high ZT re- 
quires finding a material that combines transport properties 
that do not normally occur together. Here we focus on the 
combination of high thermopower and high conductivity. 

The low T electrical conductivity of a metal or degenerate 
semiconductor depends on the electronic states and their scat- 
tering at the Fermi level, Ep, specifically a °< N(Ep) < v 2 > z, 
where N is the density of states, < v 2 > is the average Fermi 
velocity for the current direction, and X is an inverse scattering 
i ^ , ' rate. JH 01 At finite temperature the expressions are similar 
but are integrated with the derivative of the Fermi function, 
which amounts to a temperature broadening. The conductiv- 
ity therefore improves as one moves Ep away from the band 
edge, as in that case both the velocity and N(Ep) increase. The 
thermopower is different. At low T, S(T) T(da/dE)/a, i.e. 
S/T is large near the band edge where the logarithmic deriva- 
tive of a with energy is high. 

Hicks and Dresselhaus suggested that this conundrum 
could be overcome in quantum well structures. B5D They ob- 
served that in a two dimensional system the dependence of the 
density of states on energy for a parabolic band is a step func- 
tion, meaning that for the in-plane direction one may expect a 
faster onset of the conductivity with energy and furthermore 
that S will be higher for a given carrier concentration. Viewed 
in three dimensions, the Fermi surfaces of superlattices or two 
dimensional semiconductors are in the shape of cylinders or 
pipes running along the direction of the layering rather than 
the spheres or ellipsiods of three dimensional doped semicon- 
ductors. 

However, most applications of thermoelectrics involve 
macroscopic devices that are difficult to implement with su- 
perlattices and additionally, there can be problems such as 
parasitic heat conduction in barrier layers of superlattices. 
Nonetheless, it is interesting to observe that Na x Co02, which 
is representative of the highest performance oxide thermo- 



electrics, and has high ZT at high carrier concentration, [6] 
has a very two dimensional electronic structure. |0] This ma- 
terial illustrates another problem with using 2D electronic sys- 
tems as thermoelectrics. The high electrical conductivity is re- 
alized only in the layers, not perpendicular to them, while the 
heat conduction is more isotropic. As such, the very high ZT 
is realized only in single crystals for in-plane conduction or at 
least in highly textured ceramic. Here we propose an alternate 
resolution of the conundrum of high a and high S using low 
dimensional electronic structures. 

We observe that it is possible to have an electronic struc- 
ture that is low dimensional in a material that is not low di- 
mensional provided that symmetry is obeyed. This is known 
in metallic materials, the best example being body centered 
cubic Cr metal, where flat (i.e. ID) parts of the Fermi surface 
give rise to an nesting induced spin density wave. [81] An- 
other example is the superconductor Si"2Ru04, which dispite 
its tetragonal symmetry has flat one dimensional sheets of 
Fermi surface that give rise to nesting induced peaks in its 
susceptibility. Ol Generally, these cases are large Fermi 
surface metals, which are not of interest as thermoelectrics. 
However, there is no symmetry or other fundamental reason 
that this must always be the case and we begin by pointing 
out counterexamples. 

The face centered cubic rocksalt structure chalcogenides, 
PbTe, PbSe, PbS and SnTe are the basis of excellent ther- 
moelectric materials. HI 01 While the thermoelectric proper- 
ties of these materials has been discussed in terms of various 
physical models, band structure calculations in combination 
with standard Boltzmann transport theory can reproduce and 
predict their thermopowers, as illustrated by predictions for 
PbSe. 1112, 13] As is well known, the valence band (p-type) 
electronic structure is dominated by L-point hole pockets for 
low carrier concentrations and T, while at higher carrier con- 
centrations and T transport and other data imply additional 
electronic features, often discussed as a second heavy band. 
Ill4l - ll6ll According to band structure calculations, there is no 
second heavy band, but instead connections develop between 
the L-point pockets near but not at the valence band maximum. 

This is illustrated in Fig. Q] which shows energy isosur- 
faces for the near valence band edge of PbTe, PbSe, PbS and 
SnTe. These are based on calculations done with the aug- 
mented planewave plus local orbital method, II 1 711 as imple- 
mented in the WIEN2k code. [18] The calculations included 
spin-orbit, which is needed for these materials. We employed 
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We suggest that this may be a useful paradigm in the search 
for new high performance thermoelectric materials. 



FIG. 1. (color online) Calculated valence band constant energy sur- 
faces of PbTe, PbSe, PbS and SnTe, at 0.25 eV, 0.49 eV, 0.61 eV, 
and 0.41 eV below the valence band maximum, respectively. The 
corresponding carrier concentrations in holes per unit cell are 0.005, 
0.030, 0.054 and 0.016, respectively. 



the modified Becke- Johnson potential of Tran and Blaha (TB- 
mBJ), [19] which generally gives improved band gaps for sim- 
ple semiconductors and insulators. B19H21I1 Besides these de- 
tails the calculations are similar to those presented previously. 
1 12, 22- 3^1 The densities of states (not shown) show low val- 
ues characteristic of a light band up to the energy where the 
L-point pockets connect, at which point there is a sharp onset 
of a steeply rising density of states, which is clearly beneficial 
for obtaining enhanced S(T) at doping levels near the onset 
and was discussed in the context of the thermoelectric perfor- 
mance of PbTe. i27ll Here we associate this with the pipes. 

Qualitatively, the Fermi surface of a doped superlattice or 
other 2D semiconductor is cylindrical running along the stack- 
ing direction. The conductivity is low along the cylinder and 
high in the plane. Considering for example the conductivity 
along x for a cubic network of pipes running along k x , k y and 
k z as is approximately the case in these materials, the pipes 
along k y and k z will contribute as in a superlattice material in 
plane, while the pipes along k x will behave like the stacking 
direction and will not contribute to the conductivity. Thus the 
energy dependence and other behavior will be the same as the 
superlattice, including the enhanced 2D behavior of the ther- 
mopower, except that now the properties will be isotropic due 
to the cubic symmetry and superposition of pipes on different 
directions. 

Clearly, the electronic structures of the chalcogenides 
shown in Fig. Q]are approximations of this idealized behav- 
ior. Nonetheless, they suggest elucidation of the behavior of 
a cubic or other three dimensional semiconductor with a low 
dimensional electronic structure in the sense discussed above. 



CALCULATIONS 



To describe the behavior of the transport in the aforemen- 
tioned "pipes" scenario, we pursue some calculations to eluci- 
date the relevant physics. We consider the case of a one band 
material with a pipe-like electronic structure. We begin with 
the assumption that the electronic scattering time x(E) is in- 
dependent of energy, i.e. the "constant scattering time approx- 
imation" (CSTA). The CSTA has been used with quantitative 
accuracy to describe the thermo pow er of a substantial number 
of thermoelectric materials 1I271 13 11 - 43811 so that its usage is on 
solid practical grounds. 



Then we have the canonical expressions for the electrical 
conductivity o(T) and Seebeck coefficient S(T): 



a{E) = N{E)v 2 (E)x(E) 
g{T) = -J dEa(E)df(E-fi)/dE 
k B 



(1) 

(2) 



S(T) = - 



eo(T) 



dEa(E)^—^df(E-^)/dE (3) 



where / is the Fermi function, e the electronic charge, ks is 
Boltzmann's constant, i(E) is the scattering time, v(E) the 
Fermi velocity, /j the chemical potential and N(E) the den- 
sity of states. The tensor indices are suppressed for clarity, 
and in addition the integrations in actual calculations involve 
a Brillouin-zone sum. 

We now compare the thermopower and power factor S 2 <3 of 
two specific idealized Fermi surface topologies: a two dimen- 
sional cylindrical Fermi surface connecting the L-points of the 
fee Brillouin zone, as suggested by Figure 1 , and a three di- 
mensional spherical Fermi surface. This is an idealization, 
as in an actual material Fermi surfaces which contact Bril- 
louin zone faces must do so at perpendicular angles, so that 
the pipes must reconnect at the L-point pockets, as they do in, 
for example, band structure calculations for PbTe. Both bands 
are assumed to be parabolic, and to ensure a fair comparison 
we choose the radial mass of the cylinder and of the sphere 
to be equal. Additionally, as was noted by Ref. 39, in the 
chalcogenides the cylindrical band is twelve-fold degenerate 
and we have assumed this here. For comparison purposes we 
have assumed the spherical Fermi surface to also be twelve- 
fold degenerate. 

Then within the CSTA the above integrals are easily eval- 
uated for both the 3D and 2D cases, yielding the following 
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expressions (here r| = p/T, the reduced chemical potential.) 

(4) 



S3D( T ) = o F ZTx 



3^i/ 2 (Tl) 



<*3D(T) 



pe 2 l 



&D (D = 2§fi!i-T, 



02D(r) = 



^o(ri) 
2pe 2 X 
3m* 



(5) 
(6) 

(7) 
(8) 



where Af(E) is the density of states, m* the carrier effective 
mass, and F is the Fermi-Dirac integral, defined as 



Here p is the carrier density given as 

p= f dEN(E)f(E — p) 



Jo 



■Tl) + 1) 



(9) 



The 2/3 factor for the two-dimensional conductivity arises be- 
cause, each of the cylinders contributing to the density-of- 
states, and hence p, conducts in only two of three directions. 
The final piece of information necessary to calculate the ther- 
mopower is the relation of the reduced chemical potential r| 
to the carrier concentration p, which is done (for each of the 
two cases) by inverting Eq. 8; details of this procedure can be 
found in Ref. 0. 

With these mathematical preliminaries completed, we now 
move to the calculated results. For concreteness (although 
the results do not sensitively depend on these assumptions) 
we have assumed an fee cell of lattice constant 6.46 A, band 
masses of 0.2 mo, where mo is the free electron mass, and 
fixed the temperature at 1000 K. This is the approximate max- 
imum operating temperature of the chalcogenides. For the 
electrical conductivity we have assumed a doping independent 
scattering time x of 10~ 15 sec, which yields high temperature 
conductivities of 100 - 1000 (£2-cm) , in line with experi- 
mental results on these materials. In Figure 2 we present the 
calculated thermopower results for the two scenarios. As is 
evident, the 2D thermopower exceeds the 3D values by a sub- 
stantial margin throughout the entire range from 0.001 - 0.5 
holes/unit cell. At the heavy dopings of 0.05 - 0.1 per unit 
cell, the 2D thermopower is nearly double the 3D value, which 
is highly favorable for thermoelectric performance, and this 
thermopower increase comes at a conductivity reduction (Eqs. 
5 and 7), relative to the 3D case, of only one third. Given this, 
it is not surprising that the calculated 2D power factor (Figure 
3) exceeds that of the 3D case across the entire range of con- 
centration modeled, and its maximum value is two and a half 
times the corresponding 3D maximum. It is highly likely that 
correspondingly higher 2D performance (i.e. ZT), relative to 
the 3D case, would also occur. We emphasize that the notation 
2D and 3D is to distinguish the cases, but that in both cases 
we are referring to the bulk, macroscopic measurable values 
for the cubic crystal. 
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FIG. 2. (color online) The calculated thermopower for the 2D (blue 
solid line) and 3D (red dashed line) cases. Carrier concentrations are 
given in carriers per unit cell. 
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FIG. 3. (color online) The calculated power factor S 2 o in mW/m-K 2 
for the 2D (blue solid line) and 3D (red dashed line) cases. 



ANALYSIS OF ENHANCED SEEBECK COEFFICIENT IN 
2D CASE 

The results of the previous section strongly suggest that the 
two dimensional "pipe" topology is favorable for thermoelec- 
tric performance, particularly for the Seebeck coefficient, an 
indispensable ingredient of good thermoelectric performance. 
Here we provide analytic understanding of this result. 

To a first approximation, the enhanced behavior of the 2D 
system modeled here can be traced to the relatively larger 
Fermi surface volume (or, equivalently, carrier concentration) 
of a 2D cylinder relative to a 3D sphere, for given Fermi en- 
ergy. In this case the Fermi surface volume of the cylinder 
is proportional to the length of the cylinder = 2 j, which is 
a value much larger than the radius of the cylinder or the 
sphere, so that for given Fermi energy the carrier concentra- 
tion is much larger. The Fermi energy is relevant because of 
the well-known Mott formula for the thermopower, 



n 2 k B 
3e 



k B Td\og{a(E))/dE lE=EF 



and for a parabolic 3D band yields 
n 2 k B 



S = 



2e 



k B T/E F 



(10) 



(11) 
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so that the thermopower is inversely proportional to the Fermi 
energy. In two dimensions, at fixed carrier concentration the 
Fermi energy is much smaller than in three dimensions, and 
the thermopower is enhanced as a result. 

To gain additional insight into this phenomenon, as well as 
explore the effect of changing parameters such as the effective 
mass and temperature, we now pursue analytic calculations 
within two well-known limits for which closed form results 
are available: the degenerate limit, when r| = Ep /T ^> 1, and 
the non-degenerate limit, when r| < 0. High thermoelectric 
performance is typically found somewhat between these two 
regimes, but together these regimes account for most of the 
behavior of the thermopower in Figure 1 . We begin with the 
degenerate limit. In two dimensions, for radial mass m*, using 
the Mott formula it is easy to show (assuming a band degen- 
eracy of 24, 2 for spin and 12 for the 12 "pipes") that the 
thermopower takes the form 



% 2 k B 3m*a 2 k B 
3 e nh 2 p 



(12) 



where p is the carrier concentration per unit cell, and similarly 
for 3D (assuming the same band degeneracy) 



Sid 



% 2 k B 2m*a 2 k B T 



h z (% 2 p)7 

so that one finds the simple result that 

S 2 d/Sw = (ii/p)? 



(13) 



(14) 



Since the carrier concentration p per unit cell is typically 
much less than unity, one finds S2D to be substantially larger 
than S3D. Note also that for large unit cells, one requires pro- 
portionately greater carrier concentrations per cell to keep the 
same carrier concentration per cubic centimeter, so that for 
larger cells the above equation has a smaller range of validity. 
Numerically for p=0.5/u.c. (yielding an r|2D of 5.5) this ratio 
is 1.845, while the exact result is 1.747, less than a 6 percent 
difference. 

To treat the non-degenerate limit, we note that in this limit 
thermopower is generally logarithmic in carrier concentration 
(the straight lines in Figure 1), so it is the prefactors that are at 
issue. The non-degenerate limit is typically specified by r| <C 
0, so that the Fermi function (exp(^ i! ) + l) -1 reduces simply 
to exp ^y- and the integrals involving the Fermi function can 
be done exactly. In this limit, as is well known 14511 . the three 
dimensional parabolic band thermopower is given by 



e 2 



(15) 



It is a simple matter to work out the corresponding value for 
our 2D cylindrical parabolic band and one finds that 



%,r)2D = — (2-ll2D(p,r)) 



(16) 



Note that r|2D and t\t,d vary due to the topology difference, 
and we now work out an expression for their difference. For 



2 dimensions, the relation of r| and p can be evaluated easily 
and is simply (for all temperatures) 



T|2d =log(exp( : 



Tip 



1) 



1 3m*Ta 2 ' 

and in the non-degenerate limit this becomes simply 

1 f n P ■> 
ll2D=l0g( 3^ ) 



(17) 



(18) 



One can similarly work out an expression for T|3d in the non- 
degenerate limit and one finds 



r\ iD = log 



4n 3/2 r 



3(2m*)ia 3 ri 



(19) 



so that, restoring the appropriate powers of h and k^ one finds 
that 



1120-1130 = -lof 



2nh 



(20) 



For the modeled situation (to* = 0.2too,T = lOQQK,a = 
6.46A) the difference is -2.097 so that in the non-degenerate 
limit one finds S2D — Sid = l-597k B /e = 137 pV /K, which is 
very close to the difference in these values at the left hand of 
Figure 1 . This is a substantial increase, needless to say. 

The last equation reveals that if the effective mass (which 
was chosen on the basis of effective masses in the chalco- 
genides and Bi2Te3) is larger, the effective benefit in the non- 
degenerate limit is smaller, but for large effective mass materi- 
als one is typically closer to the degenerate limit. Conversely, 
if the temperature is smaller (such as for room temperature ap- 
plications) the difference is correspondingly greater, provided 
the sample remains in the non-degenerate limit. 



SUMMARY AND CONCLUSIONS 

To summarize, we have here shown that (1) low dimen- 
sional electronic structures can occur even in cubic semi- 
conductors, and that (2) such electronic structures are highly 
beneficial for thermoelectric performance. This represents 
a new paradigm for high thermoelectric performance: low- 
dimensional electronic structures enhancing performance in 
fully three dimensional bulk thermoelectrics. We suggest 
searching for new thermoelectric materials among such com- 
pounds. One such compound may be SnTe 
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